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ABSTRACT 


I.I. The PHOENIX code 


We present the results of a set of model atmospheres and 
synthetic spectra computed with the PHOENIX code. 
The models cover a range of effective temperatures 
(2TOOK < Teff < 10000K), gravities (-0.5 < log(5) 
< 5.5) and metalicities (—4.0 < [Z/H] < +0.5). In 
addition, variations of alpha elements are considered 
for each metalicity. The models are computed with a 
homogeneous set of input data in order to allow for 
direct relative comparison between the models. Eor 
example, all models use a mixing length of 1/Hp = 2.0. 
We provide synthetic spectra with a resolution of 0.2 nm 
from the UV to the infrared for all models. We give a 
brief overview of the input physics and show illustra¬ 
tive results. All synthetic spectra are available via ftp; 
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ftp.hs.uni-hamburg.de/pub/outgoing/phoenix/GAIA 


Key words: Stars; atmospheres; Stars: late-type; Radia¬ 
tive transfer; Gaia. 


I. INTRODUCTION 


Model atmospheres and synthetic spectra are of fun¬ 
damental importance in the development of filter and 
resolution systems and also helpful for the design of 
the analysis packages for the GAIA mission. Here, 
we present our first results of a grid of model at¬ 
mospheres and synthetic spectra computed with th e 
PHOENIX vl3 code (iHaiischildt & RaronI Il999h . 

The grid covers the parameter range of interest to 
the GAIA community (see Table [0, ~ 44000 mod¬ 
els and spectra in total. It is publicly available via ftp; 
ftp.hs.uni-hamburg.de/pub/outgoing/phoenix 
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parameter 

range 


step 

TeS 

2700. 

. .5000K 

100 K 


5000. 

.. 10000K 

200 K 

logg 

-0.5 .. 

. 5.5 

0.5 

[Z/H] 

+0.5 .. 

. -3.5 

0.5 

[a/a©] 

-0.2.. 

. +0.8 

0.2 


Table 1. Grid Parameters. 


Figure 1. The effect of temperature variation on a model 
with different metalicities and parameters as indicated 
in the figure. At lower temperatures (< 4000K) met- 
/dM^bii'ome important opacity sources especially in the 
UV, shifting the flux maximum of the spectrum to longer 
wavelengths, while at higher temperatures (> 7000K) 
metals are ionized such that hydrogen is dominating the 
spectrum. 

PHOENIX is a multi-purpose stellar atmosphere code 
used to compute A-T dwarfs, giants, O-B stars 
and stellar winds, exoplanet atmospheres (also irradi¬ 
ated), nova atmospheres, SNe I & II and CVs. It 
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was develoDed bv our group to be a very s 
NLTE stellar atmosphere code (iHauschildtll 19921 

sneral 

1993 

19951 lAllard & Hauschildt 1995t Ifiaron et all 

1996 

Hauschildt et al.ll 199(4 1 19971 Baron & Hauschildtl 

1998 

Allard et al.l200lHHauschildt et al,|2()C)l|) . The code can 


handle extremely large model atoms as well as line blan¬ 
keting by many millions of atomic and molecular lines. 
The radiative transfer in PHOENIX is solved in spheri¬ 
cal geometry and includes the effects of special relativity 
(including advection and aberration) in the models. The 
PHOENIX code also allows us to include a large number 
of NLTE and LTE background lines and to solve the ra¬ 
diative transfer equation for each of them without using 
simplifying approximations. Therefore, the line profiles 
must be resolved in the co-moving (Lagrangian) frame, 
which requires many wavelength points (typically 50 000 
in the GAIA grid). Details of the numerical methods used 
in PHOENIX can be found in the literature cited above. 
In order to take advantage of the enormous computing 
power and vast memory sizes of modern parallel comput¬ 
ers, we have developed a parallel version of PHOENIX 
which is now the default production version for all our 
model grids. 


2. MODEL ASSUMPTIONS 


The models of the GAIA grid are calculated for ID 
spherical symmetry and a fixed stellar mass of 1 Mq 
for simplicity. Our basic assumptions for the grid pre¬ 
sented here are time independence, hydrostatic and ther¬ 
mal equilibrium, and LTE. Convection is treated with a 
mixing length of I/Hp = 2.0. 

The micro physics setu p is similar to the one used in 
our latest model grids (lAllard et which give 

the currently best fits to the observed spectra of M, L, 
and T dwarfs. Water- an d TiO-lines are taken from the 
AM ES calcula ti on fro m iPartridge & Schwenkd Jl997h 
and ISchwenl^ ill 9981) . Our combined molecular line 
lists include about 700 million molecular lines. Typically 
about 15-300 million of these are selected in a model. 
Atomic and molecular line opacities are selected dynam¬ 
ically for the relevant LTE background lines from the 
master line lists and sum the contribution of every line 
within a search window to compute the total line opac¬ 
ity at arbitrary wavelength points. We do not use pre¬ 
computed opacity tables but individual line profiles. This 
approach allows detailed and depth dependent line pro¬ 
files to be used during the iterations. This is important 
i.e. for M dwarf spectra where line broadening and blan¬ 
keting are crucial for the model structure and therefore 
the computation of the synthetic spectra. Eor details see 
iHauschildt et alJlf2003i) . 

In the present models we have included a constant statis¬ 
tical velocity field, x = 2kms“^, which is treated like 
a microturbulence. Lines are chosen if they are stronger 
than a threshold, T = /Kc = 10“"^, where and Kc are 

the extinction coefficient of the line at its center and the 
local 5 — f absorption coefficient, respec tively. Details of 
the line selection process are given in IHauschildt et alJ 

The equation of state (EOS) is an enlarged and enhance d 
version of the EOS used in lAllard & Hauschildd il 199.51) . 
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Figure 2. By varying the gravity, we observe different ef¬ 
fects. At high temperatures (top panel) the Balmer jump 
gets more pronounced with lower gravities. At high met- 
alicities it can have an important effect on molecular 
lines as shown in the lower panel for a [Z/ H] =+0.5, 

= 3 000 K model, because the high pressure influences the 
EOS, resulting in larger concentration of molecules. The 
middle panel shows an intermediate state. 


About 1000 species (atoms, ions, molecules) are in¬ 
cluded. Eor cooler models, the formation of dust particles 
has to be considered in the EOS. In our cooler models we 
have used the assumption, that dust forms and immedi¬ 
ately rains out completely below the photosphere. Eor 
effective temperatures above 2 700 K, dust formation can 
be neglected and therefore this assumption is uncritical 
for the GAIA grid. 
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Figure 3. The variation of the spectrum with changing 
a-element abundance (compared to the solar value). The 
top panel shows a model with [Z/H] = —2.5. To em¬ 
phasize the effect of a-element abundances we show the 
difference spectra in the lower part of each panel. The 
middle panel shows a model with the same parameters 
as above but a higher metalicity [Z/H] = +0.5, and the 
lower panel a 3 000 K model. This shows that a-elements 
become more important in cool and/or metal rich atmo¬ 
spheres. Note that cool models are influenced over larger 
wavelengths ranges. 


demonstrate the behavior of the stellar spectrum. A re- 
vievy about the quality of the synthetic spectra is given in 
iKucinskas et al.l (l2005b. 

FiglU shows the effect of temperature variation on a 
model atmosphere. At lower temperatures metals become 
important opacity sources especially in the UV, shifting 
the flux maximum of the spectrum to longer wavelengths. 
At higher temperatures metals are ionized and the spec¬ 
trum is dominated by hydrogen. This effect is indepen¬ 
dent from the metalicity (upper and lower part of Fig^. 

Increasing the metalicity as shown in Fig. 0] resulting 
metal lines dominate the spectrum especially in the UV 
region. 

By varying the gravity, we observe different effects 
(Fig. m Hot, metal poor models show a more pro¬ 
nounced Balmer jump with lower gravities (top panel), 
whereas in metal rich models (middle and lower panel) 
the gravity can have important effects on the formation 
of molecular lines. The high pressure influences the EOS, 
resulting in a larger concentration of molecules. 

The GAIA grid also includes different a-element abun¬ 
dances for each set of parameters. Some are shown in 
Fig. 0 We have varied the abundance of the a-elements 
O, Ne, Mg, Si, S, Ca, Ti in our models. The effect of 
a-element abundances is small, so each panel shows also 
the difference spectra as well. Comparing the middle and 
lower panel of Figl^one sees, that a-element abundances 
become important at low temperatures and/or high met- 
alicities. 
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Figure 4. Models with Teff = hSOOA' and different met¬ 
alicity. With increasing metalicity more and more metal 
lines appear. Especially in the UV the spectrum is domi¬ 
nated by metal opacity. 


4 . CONCLUSIONS & OUTLOOK 


3. THE GAIA GRID 


In the following, we will show some examples from the 
GAIA grid. Varying one model parameter while keep¬ 
ing the others fixed, we generated a series of plots to 


The present GAIA grid provides a very detailed set of 
synthetic spectra to support the photometric working 
group. We plan a full extension of the grid to the temper¬ 
ature range between 10 000 K and 50 000 K. These calcu¬ 
lations will be done in full NLTE. It is also possible to 
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extend the present synthetic spectra to higher resolution 
upon request to the authors. 
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